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I. INTRODUCTION 



Hydrogen is the most abundant element in the Universe. It constitutes 75% of the total 
matter by mass and over 90% by number of atoms. In dense clouds of the interstellar medium 
(ISM), hydrogen exists predominantly in its molecular form and is the chief constituent. The 
hydrogen molecule is of fundamental importance in the evolution of the Universe for two 
reasons, (a) Because of its high efficiency as a coolant it increases the rate of collapse of 
low mass interstellar clouds, (b) Once ionized by UV photons or by cosmic rays, it is of 
paramount importance in all the reaction schemes that form most of the molecular species 
in the gas phase.- Moreover the hydrogen molecule enhances the sticking probability of 
hydrogen atoms and molecules on the icy interstellar dust grains.- 1 ^ 

We know by now that gas-grain reaction are the most efficient route of molecular hydrogen 
formation in the ISM.-i 5 - Two main mechanisms contribute to this formation route: (1) the 
Eley-Rideal (ER) mechanism in which a gas-phase hydrogen atom reacts directly with an 
adsorbed atom on the grain, and (2) the Langmuir-Hinshelwood (LH) mechanism in which 
both atoms are adsorbed and at least one of them diffuses on the surface of the grain to find 
its partner and to form a molecule.- Then the sticking of atomic and molecular hydrogen 
onto the interstellar dust grains plays a major role in understanding the surface chemistry 
in the ISM. 

In the dense interstellar medium dust grains (carbonaceous or silicates) are covered 
with icy mantles mainly composed of water. It is widely accepted that these man- 
tles have an amorphous structure,- - - and are covered with molecular hydrogen (1.5 - 
2.0 10 14 molecules.cm- 2 ).-^^ 1 - 

The sticking of atomic and molecular hydrogen on amorphous ice surfaces under interstel- 
lar conditions has been extensively studied theoretically over the years, but few experimental 
works have been carried out on the subject so far. The first theoretical studies were those of 
Hollenbach & Salpeter,— followed by Burke & Hollenbach,— and Leitch-Devlin & Williams.— 
But the first intensive calculations were those of Buch & Zhang. 15 These authors studied the 
sticking of H and D atoms on an amorphous water ice cluster made of 115 H2O molecules 
by using classical molecular dynamics (CMD) simulations. They calculated the sticking co- 
efficients S(E) for several beam kinetic energies E ranging from 50 K to 600 K. The fitting 
of their results was done using the simple exponential decay function S(E)=e~ E ^ E ° , where 
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E is the kinetic energy of the incident gas and where _E'o=200 K for D atoms and £'o=102 K 
for H atoms (~2 times lower than that of D). 

Matsuda et al.,— and Takahashi et a/., 17 also studied the sticking of atomic hydrogen 
on amorphous ice slab made of 1000 H 2 molecules (40Ax40Ax20A) at 10 K using CMD 
calculations. They calculated the sticking probability of hydrogen atoms as a function of the 
incident beam temperature (i.e., the temperature of the gas in the beam) and found that, 
for a kinetic energy ^=10 K, the sticking probability of a hydrogen atom is unity. 

In the same context, classical trajectory (CT) calculations were performed by Al-Halabi 
& van Dishoeck. 18 Their results on the sticking probability of H atoms on amorphous solid 
water (ASW) ice (6 bilayers of 360 H2O molecules) at 10 K were fitted by the same decaying 
exponential function used by Buch & Zhang.— These authors found that S(E) = ae~ E ^ E °, 
where E is the kinetic energy of the incident atoms, £'o=300 K (as compared with 102 K 
for Buch) and a = 1 is a constant parameter. With these fitting function and parameters, 
the sticking coefficient of an H atom with E=10 K on a surface at 10 K is equal to 0.97 and 
with £=300 K it is equal to 0.37. 

To date, only one set of experiments on the sticking of atomic hydrogen on amorphous 
water ice surfaces can be found in the literature,- and few experimental works have been 
conducted to measure the sticking coefficients of hydrogen and deuterium molecules on the 
same ice surfaces at low temperature.- 1 ^ 

Schutte et al., 3 and Govers et al.,- presented results on the variation of the sticking 
and accommodation of molecular hydrogen (deuterium) with surface coverage by H2 (D2). 
They both used bolometer experiments to study the sticking coefficient on a surface of a 
cryodeposit of H2O, N2 and Ar in the 3.5 K-15.5 K range. They found that at a surface 
temperature <10 K, the sticking coefficient of impinging H2 molecules on an initially H2 free 
surface increases very slightly, when the surface temperature decreases. This means that 
the sticking coefficient is rather independent of the surface temperature in the 3.7 K-10 K 
range. These authors also found that the sticking probability is highly dependent on the H2 
coverage of the surface, and that it increases with the increasing amount of adsorbed H 2 . 
Govers et al.fi found that the sticking coefficient of H2 is equal to 0.08 ±0.05 and that of D2 
equal to 0.27 ±0.05 for beams at room temperature and with an incidence angle of 45° with 
respect to the normal to the surface. [These values of the sticking coefficients are obtained 
directly from their figures 5 and 6. They are not exactly those indicated in the abstract of 
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their article.] 

Hornekaer et al,— studied the HD formation efficiency on both porous (p-) and non- 
porous (np-) ASW ice at 8 K using the TPD (Temperature Programmed Desorption) 
technique. These authors measured the sticking coefficient of D 2 with the King Sz Wells 
method,— shortly described in section II. They found ,Sd 2 =0.20±0.15 for a D 2 beam direc- 
tion perpendicular to the surface, at room temperature, on an np-ASW ice at 8 K. 

Amiaud et al., 21 used the same method also for D 2 and found S , £ )2 =0.38±0.05 for a D 2 
beam at room temperature and with an incidence angle of 43° with respect to the normal. 

In the present article, using the King & Wells method, we give the first experimental 
results on the variation of the sticking coefficient of molecular hydrogen and deuterium on 
ASW ice surfaces as a function of the incident molecular effusive beam temperature. 

Moreover, in order to explain: (a) the isotopic effect and the beam temperature behavior 
of our data, (b) the discrepancy between all known experimental results (including our data) 
on the sticking coefficients of H 2 and D 2 , a simple physical model involving few parameters 
is developed. To our knowledge, it is the first time that a model reproducing a so wide 
range of results is proposed. This model is divided in two parts. The first part is devoted 
to the description of the physical process of sticking alone, while the second part describes 
the effect of the effusive molecular beam velocity distribution on the measurements. 

The article is organized as follows. In section II we briefly describe the experimental 
setup and procedures. In section III we present our experimental results. In section IV 
we develop our model. In section V we apply our model to analyze our data and those of 
previous experiments (also based on molecular beams). We then discuss the obtained results 
(Sec. V) before concluding (Sec. VI). Finally the strength of our model is tested (Appendix 
A) by fitting the theoretical results of V. Buch et al— on atomic hydrogen and deuterium. 

II. EXPERIMENTAL 
A. Apparatus 

The experiments are performed using our FORMOLISM (FORmation of MOLecules in 
the InterStellar Medium) set-up, that we will briefly describe here (more details can be 
found in Amiaud et al— FORMOLISM is an apparatus devoted to study the reaction and 
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FIG. 1. Experimental setup 

interaction of various atomic and molecular beams on dust grain analogs under interstellar 
conditions. 

This apparatus (figure [Q is composed of an ultra-high vacuum (UHV) chamber with a 
base pressure of ~0.5 1CT 11 mbar, of an oxygen-free high- conductivity copper sample holder 
(1 cm in diameter), in good thermal contact with a cold finger of a He closed-cycle cryostat 
(ARS-Displex DE-204S), that can be cooled down to 8 K, and of a translatable quadrupole 
mass spectrometer (QMS Hiden HAL-3F) that is usually used in two positions, in front of 
or above the sample holder. 

The temperature is measured using a calibrated silicon diode clamped on the sample 
holder and connected to a temperature controller (Lakeshore 340). The temperature can be 
controlled to ±0.2 K with an accuracy of ±1 K in the range of 8 — 400 K. A micro-channel 
array doser (1 cm in diameter) is moved in front of the sample to expose it to H 2 vapor 
and to grow ASW ice films. The procedure for np-ASW ice film growth is already described 
under procedures (sample temperature). 

Hydrogen and/or deuterium reactants are introduced into the UHV chamber via two 
separated triply differentially pumped beam lines aimed at the sample holder. For this 
experiment we only use the first beam line, making a 62° angle with the normal to the 
surface, where the gas flows through an aluminum nozzle connected to a He closed-cycle 

5 



cryostat, in order to cool it down to a controlled temperature Tg before entering in the 
UHV chamber. A valve is located between the second and the third stage of the beam-line 
(separated by a 3 mm diaphragm) that is used in these experiments to create an effusive 
beam. In fact, opening the valve between the two stages allows us to estimate the amount of 
molecules that are diffusing from one chamber to the other out of the beam. At the entrance 
of the main chamber, also separated from the third stage by a 3 mm aperture, there is a 
flag that we use to intercept the beam, prohibiting the species to directly reach the surface 
of the sample holder. 

The effusive character of the beam at the end of the line (entry in the UHV chamber) is 
deduced from the Knudsen criterion that allows to distinguish effusive beams from free beams 
(supersonic) in the first chamber, but no direct measurement of the velocity distribution has 
been done. 

We call beam temperature (Tb) the temperature of the gas phase in the last part of the 
injection line (made of aluminum) : Tg is our tunable parameter. Let us notice for what fol- 
lows that an effusive beam is not a gas phase at equilibrium (in the sense of thermodynamics) 
and then its temperature is not defined. Then Tb is a parameter (having the dimension of a 
temperature) that characterizes the velocity distribution of the beam in the UHV chamber, 
but the convenient expression "beam temperature" is abusive (or misleading). 

B. Procedures 

In order to perform our experiments, we first move the micro-channel array doser 20 mm 
away in front of the sample holder maintained at 120 K. We then introduce H 2 vapor by 
maintaining the pressure in the UHV chamber at 2 10~ 9 mbar. With this method we grow 
np-ASW ice22r— . Previous calibrations showed that we need ~5 minutes to grow 100 ML 
(1.0 ML=10 15 cm -2 ) of np-ASW. Finally we cool down the surface temperature Tg to 10 K 
when the partial pressure of water monitored by the QMS is stable and close again to the 
base pressure in the UHV chamber. 

The sticking probability of reactive molecules on a surface is usually measured by the well- 
established beam reflectivity technique of King and Wells.— In this technique, the intensity 
of the molecules scattered from the surface, recorded by the QMS in the UHV chamber, is 
used as a direct measurement of the sticking coefficient. For this we procede as follows: 
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FIG. 2. D2 signal registered with the QMS during the 4 steps of each experiment (See text for 
detailed explanation). (A: Valve closed, B: Valve slightly open+flag, C: Valve widely open+flag, 
D: Beam directly aimed at the surface). The beam temperature Tb is 240 K and the sample at 
10 K. 



After the preparation of a stable and pure D 2 (H 2 ) beam, the np-ASW ice surface, held 
at 10 K, is then exposed to the D 2 (H 2 ) beam at a chosen temperature. The scheme of 
a typical experiment is shown in figure |2j It is divided into four steps during which the 
QMS is above the surface and then monitoring the indirect D 2 (H 2 ) signal reflected from 
the surface. (A) We first start monitoring the signal for 60 seconds without introducing any 
D 2 (H 2 ) into the UHV chamber. (B) Secondly we introduce D 2 (H 2 ) into the chamber by 
slightly opening the valve (between stages 2 and 3) but still blocking the beam and adding 
the flag in order to create a constant background pressure. This will allow us to estimate 
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the amount of molecules that are diffusing out of the beam. Intercepting the beam with 
the flag ensures that no direct molecules are reaching the np-ASW ice surface. (C) In the 
third step, the flag still intercepting the beam, we fully open the valve completely for 60 
additional seconds to obtain a full indirect flux of D 2 (H 2 ). (D) In the last step the flag is 
completely removed enabling the molecules to directly reach the surface. This fourth step 
may span up to 600 seconds and is decomposed in three stages. In the first stage (Dl), the 
QMS signal is linearly decreasing for ~ 100 seconds and then it starts increasing rapidly 
(D2) before it stablizes after ~ 100 seconds (D3). The same experiment is then repeated for 
several beam temperatures of D 2 (H 2 ), ranging from 28 to 350 K. 

After the first ~100 s, the signal starts to rise rapidly (D2) corresponding to a decrease 
in the sticking coefficient. This signal rise is obviously due to molecules that begin to desorb 
from the surface because their residence time becomes close to the time between two arrivals 
of impinging molecules.™ The plateau at the end of this step (D3) corresponds to a steady 
state regime where the number of sticking molecules is equal to the number of desorbing 
molecules. 

III. RESULTS 

This typical experiment described above is known as the King and Wells method^. It 
encompasses the four steps A to D. As in our experiments, the value of the D 2 signal in 
step (D3) is comparable to that of step (C) with the flag in the beam (see figure we have 
started all measurements at the fourth step. Then D 2 (H 2 ) irradiation begins when the flag 
is removed and the beam then directly aimed at the surface. Experimental results are shown 
in figure [3] for the deuterium beam. 

Figure [3] shows the D 2 normalized signal as function of the irradiation time for different 
D 2 beam temperature, starting from the fourth step (Dl) when the flag is removed. 

After the removal of the flag (step D starting here at t=0), the signal drops down dra- 
matically compared to step C. This signal lowering is explained by the sticking of the D 2 
molecules on the ice surface. One can also see that during the first stage (Dl), the behavior 
of the D 2 signal is highly dependent on the beam temperature Tb- 

A first remark concerns the initial value of the D 2 signal. One can clearly see that the 
lower the temperature is the lower the signal, then increasing with increasing Tb- This could 
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FIG. 3. Sticking of D2: normalized D2 QMS signal during step D of figure [2] for several D2 beam 
temperatures, ranging from 28 to 347 K. It shows how the sticking signal (stage Dl), of impinging 
molecules, varies with beam temperature Tb- The signal is normalized to its final steady value. 

be explained by the fact that for low Tg, the kinetic energy of the molecules is low, making 
collisions with the ice surface less elastic, thus increasing the probability for a molecule 
to stick onto the surface. As a result, the sticking coefficient of D2 increases when Tb 
decreases. A second remark is that, for high Tb, the signal starts with a noticeable linear 
decrease before it reaches the second stage D2. This almost linear decrease disappears at 
lower temperatures. This decrease has been explained by Govers et al..— These authors 
interpreted this decrease in the signal as an increase in the sticking coefficient induced by 
the presence of molecules already adsorbed on the surface. In fact, when a gas-phase D 2 
(H2) molecule impinges on a D2 (H2) molecule already adsorbed on the ice surface, the 
accomodation is greatly enhanced,- thus enhancing the sticking coefficient of the impinging 
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molecule. On the other hand, the disappearance of this linear decrease for lower Tg might 
be explained by the fact that at very low temperatures the sticking coefficient is already at 
its maximum and cannot increase any further. Then the signal rises slowly right from the 
start before reaching the rapid increase of the second stage. 

An absolute sticking coefficient can be derived from the curves in figure [31 The measured 
yield of D 2 (H 2 ) molecules Y(t) is the sum of molecules reflected by the surface Rf(t) (figure 
El part Dl) and a constant background factor B (figure [2, part B): 

Y{t) = B + R f {t). (1) 

The sticking coefficient S(t) is equal to the ratio between the non-reflected part of the signal 
and the incoming flux F: 

S(t) = ^-=|^. (2) 

In the steady state regime (figure [21 part D3), where t = oo, the incoming flux is equal 
to the reflected molecules. In this case F = Rf (oo) and then we can derive the sticking 
coefficient^: 

Sit) - Y{00) - Y{t) (3) 
b{t) ~ Y(oo)-B ■ {6) 

With this method we calculate the absolute sticking coefficients of H 2 and D 2 for different 

beam (gas) temperatures Tb as shown in figure 0] (the absolute sticking coefficient is defined 

as S(t = 0)). 

IV. MODEL 

In this section we present a statistical model that explains and fits these experimental 
results as well as some previously published experimental (from Govers et al.,— Hornekaer et 
al.,— Amiaud et al.—) and theoretical (from Buch & Zhang 1 -) results on sticking probabilities 
of H 2 , D 2 , H and D. 



A. Framework and assumptions 



In order to model the experiment, we consider an amorphous surface, flat on large scale 
(with respect to atomic scale), at a temperature T s . In the following T s is assumed to be 
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FIG. 4. Experimental sticking coefficient for H2 and D2 ranging from Tb = 28 K to 347 K. Plain 
circles are for H2 and squares for D2. The absolute uncertainties are equal to ±0.06 and are 
calculated by measuring the dispersion of the same sticking coefficient measured several times. 

fixed at the temperature of our experiment (T# ~ 10 K) and we omit this parameter in all 
functions to simplify notations. 

We also have a gas phase (in the half-space above the surface) composed of particles of 
mass m coming from an effusive beam that irradiates the surface. The gas phase velocity 
distribution is assumed to be that of an effusive beam (and not that of a thermal equilibrium): 
this distribution and some angle- dependent effects explain the difference between our results 
and already published experimental results. 

The heart of the model is the representation of the sticking process for which we consider 
that the sticking of the particles is essentially due to physisorption. 

First, if the surface was ordered (not amorphous) then we could assume that each imping- 
ing particle has only one possibility to stick on the surface: if it does not stick immediately 
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when it hits the surface then it is lost. Moreover in this situation of an ordered surface, 
the potential of gas-surface interaction essentially depends on the coordinate normal to the 
surface. Thus the main effect on sticking to the surface is due to the normal component v± 
of the gas particle velocity. 

In our case the surface is amorphous. This is treated as follows. Since we are interested 
in probabilities and not in the details of trajectories, we can always define the probability 
of sticking Sp for a particle (coming from the gas phase) that hits the surface (for the first 
time) at a given point P, i.e. the probability that the resulting trajectory leads to a stuck 
particle. If the surface is very disordered, the probability of sticking at the point P and 
at the point P' very close to P can be very different, removing any macroscopic physical 
meaning to Sp. So we imagine some coarse graining of the surface that divides it into a 
family of cells {C} such that each cell C can be represented by averaged local properties (for 
example an averaged normal vector nc). We define also the probability of sticking Sq for 
a particle hitting the surface into the cell C. Of course cells are sufficiently small to have 
different physical properties from place to place (in particular sticking properties) and this 
represents (statistically speaking) the disordered nature of our surface on large scale. 

Under these conditions, we can roughly transform the sticking dependence on normal 
velocity v± of the ordered situation, into a dependence of Sc on the modulus v of the 
velocity. This change is justified by the random topological nature of the surface on the 
cell scale: the macroscopic normal has no physical meaning on the cell scale and then the 
component v± is not meaningful. 

Now we assume that each cell C is characterized by a phenomenological velocity c(m, C) 
(depending a priori on the mass of the impinging particles), such that if the modulus v of 
the impinging particle velocity verifies v > c(m,C), the particle rebounds in the gas phase, 
and if v < c(m, C) the collision is sufficiently inelastic and the particle is stuck. Since our 
surface is amorphous, c(m, C) is randomly different from cell to cell . Therefore we will have 
to introduce a probabilistic distribution that represents the different values of c(m,C). 
In general the phenomenological velocity c(m, C) must not only depend on the mass of the 
impinging particle, but also on the particle-surface potential interaction. However in our 
case we are interested in the molecules H 2 and D 2 that have the same electronic properties, 
so that the particle-surface interaction must be essentially the same. 
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B. The sticking process (coefficient S(G, v)) 



S(G,v) is defined as the sticking probability of particles of given species (G) and given 
velocity. This quantity can be put in correspondence with the coefficient usually called S(E) 
(Buch & Zhang, 15 Al-Halabi & van Dishoeck— ) which is the sticking probability for a given 
kinetic energy. 

Two probabilities contribute to this sticking probability. To simplify, we consider that 
these contributions are independent and we multiply them to obtain the final sticking prob- 
ability S(G, v). 

• The first contribution Sq(G) is taken to be a characteristics of the species G (indepen- 
dent of the velocity). Focusing differently, this quantity can be roughly seen as the 
probability for a given particle of vanishing kinetic energy "put" on the surface (on a 
random place) to stick on it. 

• The second contribution P(v) is the probability to have v < c(m,C). It thus involves 
the probability distribution g(c) corresponding to the values of c(m,C). This latter is 
assumed to have the form g(c) = (f>(c/co(m))/co(m) where is a (unknown) surface- 
dependent function and co(m) is a parameter having the dimension of a velocity. The 
function must verify the condition f °° (p(x)dx = 1 (probability normalization). 

Thus P(v) writes 



Even if we do not know the function $, we can say from its definition that $ is a positive 
decreasing function such that $(0) = 1 and $(oo) = 0. 

It is noteworthy that the function $ does not depend on the molecule but only on the 
surface. The only parameters depending on the molecule being Sq(G) and co(m). 

In the case of H2 and D2, the experimental data seem to show a negligible dependence of 
c on the mass. So in what follows the mass dependence of c is neglected. 





Finally, the sticking probability S(G, v) can be written as: 




(5) 
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C. The measured sticking coefficient S(G,Tb) 



Actually, in view of the experimental results, one has to go one step further, and get 
from S(G,v) the sticking probability S(G,Tb) as a function of T B (the so-called beam 
temperature). This requires the introduction of a velocity distribution function f{v,T B ). 
The previous formula (Eq. [5]) must be changed into 

POO 

S(G,T B )= / S(G,v)f(v,T B )dv. (6) 
Jo 

Actually f(v, Tb) is not merely the Boltzmann distribution. Indeed, for a surface being 
irradiated with a direct effusive beam, a privileged direction exists and the real velocity 
distribution of the effusive beam must be taken into account. The velocity distribution of 
an effusive bearn^ f b {v,T B ) is 

f b {v,T B ) = Z~ l v b e ^bt b (7) 

where v b = v.u b > is the velocity component in the beam direction specified by the unit 
vector u b and Z is a normalization constant. In that case the gas phase is not at thermal 
equilibrium, and Tb is not the gas phase temperature, since this quantity is in fact undefined. 

If we call 9 the angle between the beam direction and the normal to the surface of our 
sample, two extreme situations can occur. In the first case 9 ~ 90° (the exact value 90° is of 
course not interesting in practice: molecules cannot hit the surface). Due to the geometry 
of the experiment, the surface can only see the central part of the effusive beam, that is the 
one- dimensional velocity distribution in the beam direction u b and thus v b = v. In the second 
case 9 = 0°. In that situation we can assume that the surface sees the full effusive beam 
velocity distribution. But we are only interested in the velocity modulus. So integration over 
spherical coordinates^ leads to the replacement v b — >• v 3 (the v 2 supplementary dependence 
arises from the volume element). Thus the angle dependent velocity modulus distribution 
f(9,v,T B ) seen by the sample writes 

f\9,v,T B ) = Z~ l v a ^e-^, (8) 

where we choose a{9) = 1 + 2cos 2 (6') to interpolate between the two extreme cases just 
discussed [because of the symmetry with respect to the sample plane, no linear term in cos 9 is 
involved ]. Z is a normalization constant obtained from the constraint J °° f{9, v, Ts)dv = 1. 
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From Eqs [5j [6] and [8], we thus obtain after some calculations 



oo 



\ v T °( m ) 

where T (m) = c^m/{2kB) and T is the Euler function. In our case 9 = 62° and a(8) ~ 1.44. 

D. Mass and temperature dependence of the sticking coefficients 

To simplify notations, we replace Tb by T, S(G,9,Tb) by Sh 2 (9,T) or Sd 2 (9,T), and 
Sq(G) by So(H 2 ) or So(D 2 ) to represent the sticking coefficients of H 2 and D 2 respectively. 
We also introduce the masses m# 2 and m D2 = 2m H2 , we then have T (H 2 ) and T (D 2 ) = 
2Tq(H 2 ) due to proportionnality to the mass. If we define the mathematical function $ as: 

1 f°° , 

$(0,x) = —p ^ u^^e- u ${^)du, (10) 



a{9)+l 
2 



we then have from Eq. [9] 



S D2 (9,T) = S (D 2 )$(9, ' 



T {D 2 ) / 



with T (-D2) = 2T (H 2 ). We can deduce finally: 



S H2 (9,T/2) = ^\S D2 (9,T). (12) 

This means that under our assumptions, the experimental data of Sd 2 (9, T) should be equiv- 
alent to that of Sh 2 (9, -|) up to a renormalization. In the remainder we call this scaling law 
a renormalization-dilation transform. 

This result is very simple and very interesting because it is independent of the form of the 
functions $ and a(9) that we use, and then it lends itself to a direct test of our theoretical 
assumptions. Then according to our model, the isotopic effect on our data is essentially 
explained by the factor 2 in Eq. f fl2|) which is just the mass ratio between D 2 and H 2 . 
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FIG. 5. Sticking coefficients: test of our model after our renormalization-dilation tranform of Eq. 
[T2l Closed circles: H2 data, open squares: D2 data transformed with Eq. [T2j These results are 
obtained with So(D2) / ^0(^2) = 1-1. 

V. DATA ANALYSIS FROM OUR MODEL 

A. Our experiment 

We use Eq. f lT2|) to represent on the same graph our experimental results for H 2 and D 2 
with the ratio So(D 2 ) / S (H 2 ) = 1.1 (value obtained after optimization). The result is shown 
in figure O This figure corroborates the renormalization-dilation transform in equation f[T2"j) . 

Of course to go further, in particular to obtain a fit of these experimental results, we need 



to fix the unknown function $. In making this choice we referred to the work of refs. (Il5f ). 



(25) and (18|) which suggests $(x) = e x2 . The corresponding transformed function $ is 



a(fl) + l 

$(0,x) = (l+x) — . (13) 
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FIG. 6. Sticking coefficients: fit of the experimental data obtained for = e~ x , Sq(H2) = 0.86, 
So(D<2) = 0.9 and To (#2) = 174 K. Closed circles represent H2 and closed squares D2 experimental 
points, respectively. 

The fits (Fig. [6]) are good for T > 80 K, but the concavity of the curves (especially for 
D2) in the domain T < 80 K does not seem to be the right one. 

To improve this feature we have used the function $(x) = (1 + x 2 )e~ x2 . When plugged 
in our model with adequate changes, it provides an excellent fit of the theoretical results 
of Buch & Zhangi^ for H and D atoms (see Appendix A). This test is interesting because 
it involves only the function S(G, v) of our model (and then the function $) and not the 
angle-dependent velocity distribution. 

The corresponding function $ is 

He x) = 1 + mx (14) 

1 ' J (l+x)W 1 ' 

with 0(B) = O.5(a(0) + 3) = 2 + cos 2 (#). In our case 0(9 = 62°) = 2.22. 

Fig. [7] shows the fits of our data obtained with this new function: the fits are seen to be 
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FIG. 7. Fits of the experimental results for H2 and D2 obtained with the function <£(x) = (1 + 
x 2 ) exp(-x 2 ). S (H 2 ) = 0.76, S (D 2 ) = 0.83 and T (H 2 ) = 87 K. 

quite satisfactory in the whole T range. We assume in what follows that the corresponding 
values S (H 2 ) = 0.76, S (D 2 ) = 0.83, T (H 2 ) = 87 K and our function $ are the physical 
parameters of the sticking process. 

B. Comparison with other experimental results 

The different experimental results (based on beam experiments) are presented in table 
U where we add our results (for the same temperature). We see immediately their lack 
of compatibility if we decide to compare them directly. But the situation is changed if 
we compare them through our model, because these experiments were made with different 
values of the angle 9 that modifies the velocity distribution seen by the sample. 

The result of Amiaud et alM is expected to be well-reproduced by our angle- dependent 
distribution, because the experiment is based on exactly the same experimental set-up (ex- 
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cept the angle) and the angular interpolating function a(8) has been partly elaborated for 
this purpose. We can be optimistic for the experimental result of Hornekaer et a/., 19 cor- 
responding to a perpendicular beam (8 = 0), because in that case we can assume that the 
sample sees the full distribution of the effusive beam, and this corresponds to one of the con- 
straints imposed to our distribution. In the case of Govers et al.,— the situation is different. 
The experiment was made with an angle 9 = 45°: some specific geometrical parameters of 
this experiment (different from our apparatus) are able to play a role in the angle-dependent 
velocity distribution seen by the sample. 

The comparison is presented in table [H The values obtained confirm the above expecta- 
tions. 

We may tentatively assume that the unique reason of the model breakdown for the Govers 
et al. data is due to a bad estimate of the velocity distribution parameter a. From Eqs. [11] 
and HU we deduce that the sticking coefficient obtained by Govers et al. must verify 

where So and To are our values (for H2 and D2), but (3 is now an unknown coefficient 
(identical for H 2 and D 2 ). 

Using the results of Govers et al. for H2 and D2, we can solve Eq. [151 an d we obtain two 
values of (one for H2 and one for D2). If our assumption is correct these values must be 
very close. 

Taking into account the uncertainties on the sticking coefficients, we obtain (3(H.2) G 
[2.78,4.01] and /3(D 2 ) G [2.70,3.23]. Now recalling that the full distribution of the effusive 
beam corresponds to /3 = 3 which is a physical bound, we can restrict the previous intervals 
to /3(H 2 ) G [2.78,3] and (3(D 2 ) G [2.70,3]. Then the mean values are /3(H 2 ) = 2.89 and 
/3(D 2 ) = 2.85. Thence /3(H 2 ) — /3(D 2 ) — 2.87 is compatible with our assumption stipulating 
a common value of (3 (the value obtained with our angular distribution was /3 = 2.5). With 
this value (3 = 2.87, we obtain the sticking coefficients Sh 2 = 0.12 and Sb 2 = 0.28: those are 
now in good agreement with the data of Govers et al.. 

We conclude that these different comparisons confirm both our physical model for the 
sticking process, and the description of the beam angular effect in these experiments (though 
our angle- dependent velocity distribution is not universal). Moreover this analysis shows that 
these different data must not be directly compared, and this explains the apparent discrep- 
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0.37 


Hornekaer 


D 2 


0° 


300 


0.2 ±0.15 


0.25 


Amiaud 


D 2 


43° 


293 


0.38 ±0.05 


0.36 


This study 


H 2 


62° 


300 


0.28 ±0.06 


0.24 


This study 


D 2 


62° 


299 


0.42 ±0.06 


0.43 



TABLE I. Comparison between the sticking coefficients obtained in different experiments at room 
temperature (Govers et air Hornekaer et al.^L Amiaud et a/.,? 1 our results) and our model. 

ancy between them. This analysis implies also that these experimental sticking coefficients 
are different from the true thermal coefficients (probability of sticking with a gas phase at 
thermal equilibrium) suitable for astrophysical applications. 

VI. DISCUSSION 

We are perfectly aware that the detailed interaction (and behavior) of atoms and 
molecules with the ASW ice surface are not similar. Moreover in each situation (atomic or 
molecular), we omit a lot of parameters. However one of the main interests of this kind of 
simple statistical model is precisely to erase the details and to keep the main features. The 
model could have been more complete if we were able to deduce (or compare) the velocities 
c directly from other results. 

Choosing the dependence of sticking on velocity rather than on kinetic energy can be 
justified by the following arguments. On one hand the gas-surface interaction is based on 
short-ranged potentials that implies a finite distance L of interaction, and on the other 
hand sticking implies an energy transfer to the surface that needs a minimal time tq of 
interaction. If the impinging particle has a velocity v, it interacts with the surface during 
a time r ~ 2L/v. Then the particle sticks to the surface if r > r that is v < Cq = 2L/t . 
Then there exists a characteristic velocity Cq = 2L/r relative to sticking. 

From the fits achieved for the atomic and molecular species (Appendix A and Sec. IV Aj) . 
we can compare the velocity parameters c (H) and c (H 2 ). Since atomic and molecular 
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hydrogen have not the same electronic structure, we expect cq(H 2 ) ^ cq(H). From the values 
m H cl(H) = 2E = 104 K (Appendix Ej) and 0.5m H2 c 2 (H 2 ) = m H c 2 (H 2 ) = T (H 2 ) = 87 K 
(Sec. IV Ap we obtain 

^o.»i. da) 

As expected c Q (H 2 ) ^ c (H), but the values are close. If we use our estimate c =2 L/r 
where L is a maximal distance of interaction and tq a minimal time of interaction, and if 
we assume roughly that t (H 2 ) ~ T o(H) (we have not reached a regime of limited energy 
transfer), we find 

HH- 91 ' (i7) 

This result is compatible with the idea that the range of H 2 -surface interaction is shorter 
than that of the H-surface interaction: it is more difficult to trap H2 than to trap H in a 
potential well on the surface, because H2 has internal motions with possible energy transfers 
between the different degrees of freedom of the molecule (rotation-translation of the two 
atoms). Thus if the two atoms of the molecule approach the critical interaction distance 
L(H) where only one H atom is trapped, this does not guarantee that the molecule will 
be trapped in the potential well and might be kicked out from the surface. This situation 
cannot take place in the case of only one H atom impinging on the surface. 

The angular dependence of the velocity distribution seen by the sample, and then of the 
sticking probabilities measured (for experiments based on effusive beams), is a key point to 
understand the discrepancy between the different experimental values. As a consequence, 
these values are not the "thermal" sticking coefficient (probability of sticking with a gas 
phase at thermal equilibrium). Astrophysicists can be interested only in this coefficient. It 
can be deduced from Eqs. EJ HH and [TH by remarking that the thermal velocity modulus 
distribution (Boltzmann law) takes the form v 2 multiplied by the Boltzmann exponential 
factor. So this corresponds formally to the case a{9) = 2 of our angle-dependent velocity 
distribution (EqJHJ). We then deduce (3(9) = 2.5 from Eq. [TU The equation [TT1 gives the 
final result. Then the true thermal coefficient is given (with our notations) by 

with = 2.5. 

Another point can be analyzed from our model: it concerns the sticking probability of the 
molecule HD. It is reasonable to extend the mass law obtained for H2 and D2 to the case of 
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So Tp(K) 

H 1 52 

D 1 104 

H 2 0.76 87 

D 2 0.83 174 

HD 0.8 130.5 



TABLE II. Table of the different coefficients obtained in this article. 

HD. We deduce that the temperature parameter T (HD) verifies T (HD) = (3/2)T (H 2 ) = 
130.5K. Moreover since S {R 2 ) = 0.76 and S (D 2 ) = 0.83, we can estimate S {HD) ~ 0.8. 
With these parameters we obtain the thermal HD sticking coefficient from Eq. [TH 

The values of the parameters So and To obtained for H 2 and D 2 (from our experiments: 
Sec. IV All , for H and D (from V. Buch computations: Appendix A) and for HD (prediction) 
are summarized in the table [III 

VII. CONCLUSION 

In this paper we have reported a set of experiments that we have conducted to measure 
the sticking coefficients of molecular hydrogen and deuterium on np-ASW ice surfaces held 
at 10 K using the King & Wells method. A study of the variation of the sticking coefficients 
with the molecular beam temperature was also presented. To our knowledge, this is the 
first experimental work that measures the sticking coefficient of hydrogen and deuterium 
molecules as a function of the beam temperature. 

We have also presented an original model that explains the isotopic effect and the temper- 
ature behavior of the obtained experimental data. The model succeeds in fitting the present 
data for molecular species (H 2 ,D 2 ) as well as previous experimental results on (H 2 ,D 2 ) and 
theoretical ones of Buch & Zhang on atomic species (H,D). To our knowledge, it is the first 
time that a simple model reproducing a so wide range of results is proposed. 

Following our model, we determined the values of the physical parameters describing the 
sticking process independently of any specific experimental situation (So, T , function $). 
Then from these values, we have been able to propose an explicit formula giving the true 
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thermal sticking coefficient (Eq. [T8|) for H2, D2, H and D. 

The astrophysical interest of this article is that we can now extract some useful values of 
the thermal sticking coefficients of H2 and D2 (Eq. [181 Tab. II) relevant to the dark clouds 
of the interstellar medium. 
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Appendix A: Testing our model on atomic sticking probabilities of H and D 

In their theoretical article on the sticking probabilities of H and D atoms on clusters of 
amorphous water ice (cluster temperature T ~ 10 K similar to our surface temperature), 
Buch & Zhang 15 give a set of values for Sh(E) and Sd(E). These quantities correspond to 
our coefficient S(G,v) because E = (l/2)mv 2 . Moreover the function $ of our model only 
depends on water ice properties and then it must be roughly the same. These data allow us 
to get rid of the average effect on the sticking coefficients due to our experimental velocity 
distribution and we can test directly the expression of S(G, v) and our function $. From 
our model (Eq. |5]) we have 



Let us notice that the velocity parameter cq involved in this equation is not necessarily the 
same as the one previously defined for the molecular case, because atomic and molecular 
hydrogen have not the same electronic properties. Using the procedure of Sec. IIVDI we 
obtain the scaling law (renormalization-dilation transform) 




(Al) 



S H (E/2) 



S (H) 
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FIG. 8. Sticking of H and D from Ref. I la . Closed circle for H, closed square for D. Open triangle: 
D data points transformed with Eq. IA2I (Sq(H) / Sq(D) = 0.95). Fit of H and D data points based 
on Eq. EHJ S (H) = S (D) = 1 and E = 52 K. 

where the factor 2 is the mass ratio between D and H. The effect of Eq. IA2I on their data 
can be tested independently of any fit as in Sec. IV Al This is shown on Fig. [SJ 
Using our function = (1 + x 2 )e~ x2 , we have 

S H (E) = So(H)FH) 
S D (E) = S Q (D)F(^) ' 

where F(x) = (1 + x)e~ x and E = E (H) = (l/2)m#CQ. The corresponding fits for atomic 
hydrogen and deuterium data are given on Fig. [HJ They confirm that our choice of the 
function $ is very satisfactory. Moreover this test, made on a complete different kind of 
data, shows the strength of our model and confirms the validity of our general physical 
assumptions. 
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